NORDITA-2001/11 HE 
SPIN-2001/15 



Supersymmetric 3-branes on smooth 
ALE manifolds with flux t 

o 

\ M. Bertolini 1 , V. L. Campos 2 , G. Ferretti 2 , 

p 1 P. Fre 3 , P. Salomonson 2 , M. Trigiante 4 

<D ' 

fx. ■ 1 NORDITA, Blegdamsvej, 17 2100 Copenhagen 0, Denmark 

(N 



> 
OO 



X 



2 Institute for Theoretical Physics - Goteborg University and 
Chalmers University of Technology, 412 96 Goteborg, Sweden 



3 Dipartimento di Fisica Teorica, Universita di Torino, & INFN - Sezione di Torino 

via P. Giuria 1, 1-10125 Torino, Italy 

4 Spinoza Institute, Leuvenlaan 4, Utrecht, The Netherlands 

o 



Abstract 



We construct a new family of classical BPS solutions of type IIB supergravity describing 
3-branes transverse to a 6-dimensional space with topology M 2 x ALE. They are characterized 
by a non-trivial flux of the supergravity 2-forms through the homology 2-cycles of a generic 
smooth ALE manifold. Our solutions have two Killing spinors and thus preserve J\f = 2 
supersymmetry. They are expressed in terms of a quasi harmonic function H (the "warp 
factor"), whose properties we study in the case of the simplest ALE, namely the Eguchi- 
Hanson manifold. The equation for H is identified as an instance of the confluent Heun 
equation. We write explicit power series solutions and solve the recurrence relation for 
the coefficients, discussing also the relevant asymptotic expansions. While, as in all such 
M = 2 solutions, supergravity breaks down near the brane, the smoothing out of the vacuum 
geometry has the effect that the warp factor is regular in a region near the cycle. We interpret 
the behavior of the warp factor as describing a three-brane charge "smeared" over the cycle 
and consider the asymptotic form of the geometry in that region, showing that conformal 
invariance is broken even when the complex type IIB 3-form field strength is assumed to 
vanish. We conclude with a discussion of the basic features of the gauge theory dual. 
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1 Introduction 

After the seminal paper by Maldacena fij, many efforts have been devoted to extend the 
gauge/gravity correspondence to less supersymmetric and non-conformal cases. In this context 
considerable attention was recently directed to the study of fractional branes Q-Jl^]. These are 
the natural elementary D-branes occurring whenever string theory is reduced on a (not neces- 



sarily compact) orbifold flq] , fl7|1-|20l. Because of their nature they cannot move away from the 
orbifold apex and thus the dual gauge theory on their world-volume lacks the relevant moduli 
fields. Generically, this leads to both reduced supersymmetry and non-vanishing /^-functions. 
Most interesting are the fractional D3-branes, namely the case when the world- volume theory is 
four-dimensional. In this respect the two most appealing situations are provided by the N = 1 
case emerging from singular limits of CY spaces and the N = 2 case arising from the singular 
limit of ALE spaces. Much work was devoted to both. 

A common feature of many supergravity solutions representing non-conformal situations is 
the presence of naked singularities of repulson type pi] ]. These correspond to IR singularities 
at the gauge theory level and one expects that they should be resolved or explained by some 
stringy effect. Although a general recipe does not seem to exist, progress in understanding such 
an issue was made both for the M = 1 and the N = 2 case. 

In constructing supergravity solutions corresponding to fractional branes one proceeds in 
two steps. Firstly, one defines the vacuum solution as a space-time obtained by quotiening 
ten-dimensional Minkowski space-time modulo a discrete group with all other fields set to their 
constant background values (typically, the latter are zero with the exception of the dilaton, which 
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sets the string coupling, and the NS-NS 2-form B, which determines whether the conformal field 
theory is well defined or not |22]| ). Secondly, one places the branes at the orbifold/conifold apex, 
which has the effect of turning on the relevant gauge fields and of warping the metric in a way 
that makes the curvature non-zero even away from the singularity. As already stressed, these 
solutions generally display singularities (see || for the M = 1 case and 0] for the N = 2 case) . 

Once confronted with superstring theory on such a space-time, it is a natural question to ask 
whether this space can be thought of as a singular limit of a smooth geometry. This can always 
be achieved when no branes are present. If the smoothing out can be carried out even in the 
presence of branes, and if the warp factor is such that the curvature remains small everywhere 
and the singularities in the solution are removed, then one can claim that the supergravity 
solution provides a dual description of the gauge theory. This is what was accomplished, for 
instance, in the M = 1 deformed conifold solution of Klebanov and Strassler B, where the 
dynamics of the brane deforms the conifold at its apex, yielding a singularity free solution. 

What is somewhat frustrating is that, contrary to what one would naively expect, the N = 1 
case turned out to be more amenable to a smooth supergravity analysis than the N = 2 case. 
Actually, for the N = 2 completely different phenomenon has been shown to occur 

[JjJ-ljn]]: the singular region is excised by the enhangon mechanism [ ffiH] , which corresponds to 
a locus where extra massless degrees of freedom come into play, thus completely changing the 
geometry in the interior. From the gauge theory point of view the enhangon is the scale at which 
the gauge coupling diverges and where instantonic effects become relevant. 

Despite this specific feature of the N = 2 case, one could anyhow try to see if better behaved 
solutions can arise when considering, analogously to what has been done in the conifold case in 
H and m , the blow-up of the singular space, namely the smooth ALE manifold. Then one could 
try to see if there is any relation between the enhangon and the geometric moduli of the better 
behaved geometry, try to disentangle the singular properties of the solution, and eventually 
investigate their meaning in the context of the dual gauge theory. 

Starting from the singular space C 2 /r, where V is a discrete subgroup of SU(2), one should 
inquire whether one can soften the vacuum geometry to a smooth ALE space. It is well known 
that this cannot be done in a supersymmetric way if fractional D3-branes are present and if the 
string theory is really at the orbifold point, that is, with a non zero background l?-field threaded 
through the vanishing cycle of the orbifold. There are many arguments for such a "no-go" 
theorem. Probably the most intuitive one involves going to the T-dual type IIA picture, |2jJ]. 
There one finds two NS5-branes with a stack of D4-branes (dual to the fractional D3-branes) 
stretched between them. The blowing up of the orbifold corresponds to moving the NS5-branes 
in the direction transverse to the D4-branes, and this involves stretching the D4-branes even 
further, thus increasing the energy of the system and hence breaking supersymmetryQ. 

Yet the above argument does not rule out the existence of supersymmetric bulk solutions of 
type IIB supergravity that smooth out the quotient. Such solutions, which eventually have to 
be interpreted as suitable bound states of some D-brane state, are precisely those that we find 
in the present paper. 

In particular, we find M = 2 BPS solutions corresponding to 3-brane states where the 
transverse space is given by M. 2 x ALE and where there is a non zero value of the complex type 
IIB supergravity 3-form in the transverse directions (note that this is the distinctive feature 
of fractional D-branes in the singular orbifold theories). As already noticed, solutions of this 
type with ftf = 1 supersymmetry have appeared in the recent literature but they were mostly 

1 We thank R. Russo for a discussion on this point. 
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concerned with the case where the six-dimensional transverse space was a cone. This had the 
simplifying feature of making the warp factor depend only on one (radial) coordinate. This 
feature cannot be preserved if one imposes M = 2 supersymmetry since the required SU(2) 
holonomy forces the geometry of the transverse space to be of the factorized form described 
above. The price to pay is that the warp factor is now dependent on (at least) two variables 
and it is governed by a complicated partial differential equation. In fact the whole set of IIB 
equations can be reduced to the solution of such an equation for the warp factor, whose source 
is essentially dictated by supersymmetry up to an arbitrary analytic function j(z). 

Supergravity alone is not sufficient to determine 7(2) or the boundary conditions. This 
arbitrariness implies that our solution describes various deformation or various vacua of M = 2 
theories. The case 7(2) = const, describes a vanishing three-form flux and corresponds to the 
well known conformal J\f = 2 theory with product gauge group U(N) x U(N), hyper-multiplets 
in the bi-fundamental representation and Fayet-Iliopoulos terms describing the ALE moduli. 
The new ingredient in our construction is the following. We show how at the supergravity level 
a three-form flux can be turned on consistently with supersymmetry. Some consequences of this 
for the gauge dual theory are analyzed in the concluding section yet a complete analysis of such 
an issue involves a systematic study of possible supersymmetric world-volume actions of branes 
on arbitrary curved backgrounds (in particular ALE manifolds) which is work in progress and 
the subject of a forthcoming publication. 

The paper is organized as follows. In section [2] we briefly recall the structure of the equations 
of motion of type IIB supergravity, fixing our conventions for the fields. In section ^, after 
studying the general set-up for arbitrary ALE spaces, we give the structure of our 3-brane 
solution whose distinctive feature is the localized 3-form flux along any exceptional compact 
2-cycle of the ALE space. Finally, applying a criterion which is analogous to that discussed in 
[^]-|28|], (see also ^9|) we give the proof of its supersymmetry (in this respect, our solutions 
fall in the so called type B class f30|| , characterized by the 3-form to be imaginary self-dual). In 
order to consider an explicit example, in section || we specialize to the simplest such space, the 
Eguchi-Hanson space, and study in detail the properties of the warp factor. Finally, in section || 
we summarize our results, discuss the relevant physical properties of our solutions, discuss some 
basic aspects of their gauge theory dual and conclude with some comments on open issues. 



2 Bosonic action and field equations of type IIB supergravity 



As it is well known, type IIB supergravity does not have any conventional supersymmetric action. 
However, as it happens for all on-shell supergravity theories, the complete set of field equations 
can be obtained as consistency conditions from the closure of the supersymmetry transformation 
algebra |3l|]. In the case of type IIB supergravity, one was also able to obtain a complete, 
manifestly SU(1, l)-covariant formulation of the theory based on the rheonomic approach to 
supergravity theory |33| . 

The bosonic part of the equations can be formally obtained through variation of the following 
action 0: 



Sim 




d 10 x 



-2y/-detg R 



dip A *dip + e^Ff A *F$ S 



1 „2v3 rpRR a . pRR 
+ e ^[1] A *^[l] 



2 Note that our R is equal to —^R old , R old being the normalization of the scalar curvature usually adopted in 
General Relativity textbooks. The difference arises because in the traditional literature the Riemann tensor is 
not defined as the components of the curvature 2-form R ab rather as —2 times such components. 
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where: 



+ e ^[3] A *^[3] + 9^5] A *^[5] - G [4]A^ [3] A ^ [3] 



F [1] R = dC [0] 
F [3] S = dB [2) 
F [3] = dC m ~ C[o] dB [2 ] 



(2.1) 



F3* = dC [4] -±{B [2] AdC [2] -C [2] AdB [2] ) (2.2) 



It is important to stress though that the action Q2.1| ) is to be considered only a book keeping 

7 Rt 

[5] 



device since the 4-form Cw is not free, its field strength F3, R being subject to the on-shell 



self-duality constraint: 

'[5] = **[5] 

^From the above action the corresponding equations of motion can be obtained: 



FjS* = *F& (2.3) 



d*d<p-e^F™A*F { f = -h (e-*F$ a McFff -e*F™A*F§f) (2.4) 

d(e^*F™) = -WfAtff (2.5) 

d(e^*^)+e^A*J^ = -FffiAFff (2.6) 

d(e^F™) = -F^AFff (2.7) 

d + Fffi = -Ff AF,f (2.8) 

1 e 2 ^ 
- 2 R M n = ^d M ipd N if + —d M C[ ]d N C[ ] + 150F[ 5 ] M ....F [5] ^ 

+y[ e ^[3]M--^[3]JV + e ^[3]M-^[3]iV 

'' { r e-*F» s F& s ~ + e^F™ f£P" ) (2.1 



It is not difficult to show, upon suitable identification of the massless superstring fields, that this 
is the correct set of equations which can be consistently obtained from the manifestly SU(1, 1) 
covariant formulation of type IIB supergravity [p2|. 



3 3-brane solution with ALE flux 

In this section we provide the BPS solution corresponding to a 3-brane transverse to a smooth 
ALE space, namely we construct type IIB supergravity solutions describing 3 branes on a vacuum 
IR 1 ' 3 x IR 2 x ALE. This will be achieved without an analysis of the specific form of the world- 
volume action of the brane, i.e. of the source term. Our physical assumption will just be 
that, together with the usual RR 5- form flux, the 3-brane solution has a non-trivial flux of the 
supergravity 2-form potentials along (one of) the compact two cycle(s) of the blown-up orbifold 
(this translates into a non- vanishing value of the complex 3- form field strength). 
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3.1 Solution of the bosonic field equations 

We separate the ten coordinates of space-time into the following subsets: 

fi = 0, 1, 2, 3 coordinates of the 3-brane world volume 
x M = ^ z = x 4 + ix 5 complex coordinate of M. 2 ~ C (3-1) 
y T t = 6, 7, 8, 9 real coordinates of the ALE 4-space 

and we make the following ansatz for the metric^: 

ds 2 = (-rj^dx' l dx' / )+H^dzdz + H^ds 2 ALE (3.2) 

H = H(z,z,y) 

ds\ LE = gt LE dy T dy° 

det(g) = Hdet(g ALE ) 

i]^ = diag(+, -,-,-) (3.3) 

where g^ E is the metric of any ALE space and we denote the six-manifold spanned by z, z and 
V T by 

M 6 = M 2 x ALE (3.4) 
In terms of vielbein the ansatz ( |3,2| ) corresponds to 

A _ jV a = H~ 1 / i (z,z,y)dx a , v 

V -\V*=HV\z,z,y)e i 1 j 

where e l are the vielbein 1-forms of the manifold Mq. The structure equations of the latter areQ: 

= de l - Q ij A e k rj jk 
rv = dQ ij - Q ik A Q ej flu = R ij lm e e A e m (3.6) 

The relevant property of the Mq metric that we use in solving the Einstein's equations is that 
it is Ricci flat: 

R tn L = o (3.7) 

What we need in order to derive our solution and discuss its supersymmetry properties is the 
explicit form of the spin connection for the full 10-dimensional metric (3.2) and the corresponding 
Ricci tensor. From the torsion equation one can uniquely determine the solution: 

uj ah = 

uj ai = I H^l 2 dx a rf k d k H 



'J 



Qij + hJi . = -\H~ l e [i r^ h d k H (3.8) 



3 As explained in appendix A, the conventions for the gamma matrices and the spinors are set with a mostly 
minus metric dr 2 . In the discussion of the solution, however, we use ds 2 = —dr 2 for convenience. We hope this 
does not cause any confusion. 

4 The hats over the spin connection and the Riemann tensor denote quantities computed without the warp 
factor. 
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Inserting this result into the definition of the curvature 2-form we obtain: 



R a b 

It; 



H- 3 / 2 n H H - H' 5 ' 2 d k Hd k H 



where 





-H- 3 / 2 a 6 H5i - -H- 5 / 2 d k Hd k H5 i i + ^H'^dMH 
8 J 8 " J 4 



□r2 + Dale 

is the Laplacian on the tensor product space ( |3.4| ). The equations of motion for the 
ip and C[ ] and for the 3-form field strength ^^J S and -Fjfj^ can be better analyzed 
a complex notation. Defining 

H± = ±2e-^ 2 F^ s + i2e^ 2 F§ ] R 



eq.s (2.4)-(2.7) can be written respectively as: 



while eq. ( |2.8[) for the 5-form becomes: 



iF| R A W H 



Besides assuming the structure ( |3.2[) we also assume that the two scalar fields, 
dilaton ip and the Ramond-Ramond 0-form Cr i are constant and vanishing: 

ip = ; Cfqi = 



(3.9) 



(3.10) 

scalar fields 
introducing 



(3.11) 
(3.12) 

(3.13) 
(3.14) 

(3.15) 
namely the 

(3.16) 



As we shall see, this assumption simplifies considerably the equations of motion. 

The basic ansatz characterizing our solution and providing its interpretation as a 3-brane 
with a flux through a homology 2-cycle of the ALE space is given by the following: 



H+ = 2drii{z,z) A uj 1 
H- = -2drji{z,z) A lo 1 



(3.17) 



where 7/(2, z) is a complex field depending only on the coordinates z,z of R 2 , while co 1 (I = 
1, . . . , k) constitute a basis for the space of square integrable, anti-self-dual, harmonic forms on 
the ALE manifold. Notice that upon the ansatz ( p. 16 ), H.+ equals the usual complex 3-form 
G[3| made up by the R-R and NS-NS type IIB 3-forms. More precisely 7i + = 2 Gm while 
TL- = —2 Gj 3 j. From eq.( 3.17| ) one then sees that our solutions fall in the same class of solutions 
discussed in [30], dubbed type B solutions, which are dual to M-theory solutions found by Becker 
and Becker in [34]. The former are in fact characterized by the 3-form being a self-dual 
(2, 1) form. 



G 



As it is well known f\ a smooth ALE manifold, arising from the resolution of a C 2 /T singu- 
larity, where T C SU(2) is a discrete Kleinian group, has Hirzebruch signature: 



r = (#of conjugacy classes of T) — 1 = rank of G(T) (3.18) 
In the above formula G(r) is the simply laced Lie algebra corresponding to T in the ADE 



classification of Kleinian groups, based on the Mac Kay correspondence [35|. As a result of 



eq.(3.18) the ALE manifold that is HyperKahler admits a triplet of self-dual 2-forms that are 



non-integrable and exactly r integrable anti-self-dual harmonic 2-forms. For these latter one can 
choose a basis to 1 that is dual to the integral homology basis of 2-cycles S/ whose intersection 
matrix is the Cartan matrix C of G(T). Explicitly we can write: 





= & 


JY. k 






= -c- 


u) 1 A uj j 


= -A 



II J 

ALE 

IJ 



A 1J (y)n ALE (3.19) 

where C IJ is the Cartan matrix of the corresponding (non-extended) ADE Dynkin diagram and 
A IJ (y) is a positive definite matrix whose entries are functions of the ALE space coordinates y's. 
The anti-self-duality of the or guarantees that J lo 1 A -kuo J is positive definite. Finally, following 
the conventions of appendix [A] (see in particular eq.( [A.5| ) ) Qale is the volume 4-form of the 
ALE manifold. 

If we insert the ansatze ( |3.16 , 3TT7| ) into the scalar field equations (2.4, 2.5) we obtain: 



H+ A *H+ = (3.20) 

which in turn implies that d z ^i d^j = 0. This equation is solved by choosing jj to be holomor- 
phic: 

d _ d 

dji = where d = dz — ; d = dz — (3.21) 

oz oz 

Next we consider the self-dual 5-form which, because of its definition, must satisfy the 
following Bianchi identity: 

Our ansatz for F3P is the following: (The volume forms O are described in appendix A) 



dF,§P = ihH+ A H- (3.22) 



F { f = a{U + *U) 

U = d [H" 1 O r i,3) (3.23) 

where a is a constant to be determined later. By construction F^ R is self-dual and its equation 
of motion is trivially satisfied. What is not guaranteed is that also the Bianchi identity ( |3.22| ) is 
fulfilled. Imposing it, results into a differential equation for the function H(z,z,y). Using the 
Hodge duality relations ( |A.5 ) listed in appendix [A| we easily get: 



★ U = i dH A FLale — i dH A VLale + ^ (★ d A LEH) Adz Adz (3.24) 



5 See for instance Jl(| for an early summaryof ALE geometry in relation with superstrings and conformal field 
theories. This relation was developed in Jn], [H| and is of primary relevance in connection with D-branes. 
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Taking the exterior derivative of ( |3.24j ) we get: 

d-kU = 2i BdH A VLale + ~ {d(ALE)* d(ALE)H) Adz Adz 



(□ R2 H + U ALE H) l -dzAdzA n ALE (3.25) 



Since dF RR = ad*U and 



H.+ A H- = - 4 <9 2 7/ 3 5 7j A' J dz Adz A Vt A LE (3.26) 
we conclude that the Bianchi identity ( 3.22| ) is satisfied by our ansatz if: 



(□ R2 + Dale) H = --d zlI djfu A IJ (y) (3.27) 



a 



This is the main differential equation the entire construction of our sought for 3-brane solution 
can be reduced to. We are going to show that the parameter a is determined by Einstein's 
equation and fixed to a = 1. 

Let us consider next the field equation for the complex three-form, namely eq. fl2.6|) and 
fl2.7p . Since the two scalar fields are constant we have: 

d*H + = iF [ R ] R AH+ (3.28) 

Using our ansatz we immediately obtain: 

d * H + = -2 \H~ 2 dH A fry/ Aw' A H i, 3 + 2i H^d A A n R i, 8 
iF^f A H+ = -2mH- 2 dH A &yi Aw' A n R i, 3 (3.29) 

Hence if a = 1, as we have announced, the field equation for the fields 7/(2) reduces to: 

dA<97/ = □ R a7/ = (3.30) 

This equation has to be appropriately interpreted. It says that 7/ are harmonic functions in two- 
dimensions as the real and imaginary parts of any holomorphic function 7/(2) certainly are. The 
bulk equations do not impose any additional constraint besides this condition of holomorphicity. 
However, in presence of a boundary action for the 3 brane, the equation will be modified into: 

□ R 27/=j/(z) (3.31) 

ji(z) being a source term, typically a delta function. In this case 7/ is fixed as: 

jl(z) = [ G^z^j^dz' (3.32) 



where G R 2(z,z') is the Green function of the M? Laplacian in complex coordinates and turns 
out to be proportional to logz. From this point of view one could already expect the enhangon 
mechanism to be at work in a very similar way as in the orbifold case J7[,|15[|. We will come 
back to this issue in section 4. 

In order to show that also the Einstein's equation is satisfied by our ansatz we have to 
calculate the (trace subtracted) stress energy tensor of the five and three index field strengths. 
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For this purpose we need the components of These are easily dealt with. Relying on the 

ansatz ( UD and on eq.(|3.5|) for the vielbein we immediately get: 



Fiabcd 

F 

otherwise 



A\ fi e 



abed 



5l e *ii— Js 



where: 



fi = -aH- 5 / 4 diH 
Then by straightforward algebra we obtain: 



(3.33) 



(3.34) 



rpl 

j 



t?RR 
F [5] 



-75F a ---F b .... = ~5 a b f p f 
-a 2 l5iH-^' 2 d r> Hd p H 



RR 



[5] 



-75 F l 



Fa 



iffj ~ ~J)fvf P 



a 



4 8 J 



(3.35) 



Let us next consider the complex 3- forms TL±. With our ansatz (13.171) they satisfy the 6- 
dimensional self-duality relation on Mq: 



*6 7~t± = ii7i± 
which, in components, is rewritten as follows 



±iW± 



ijk 



1 

3! 



T-/ tr ' 



Iran 



(3.36) 
(3.37) 



Eq.( |3.37 ) implies that: 
2J[W] = 



1 



(3.38) 



where we have used the result 
H+ A 



= 3!H + ""H_|... x Vol 
having named the complete volume form of the ten-dimensional space-time: 

Vol = H 1 ' 2 O a i )3 A Oals A n R 2 = ^-detgdX A ... A dX s 
Inserting eq.s ( |3.35| ), fl3.38| ) and ( pT9| ) into Einstein's eq.(|2.9[): 



(3.39) 
(3.40) 



rpd 

1 b 



J 



F, 



RR 



F, 



[5] . 
[5] 



3 

+ T] [H] 



(3.41) 
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we see that they are satisfied, provided 

a = 1 (3.42) 



and the master equation ( 3.27|) is satisfied. This concludes our proof that an exact 3-brane 



solution with a non-trivial flux through a homology 2-cycle of a generic smooth ALE manifold, 
does indeed exist. In the next subsection we prove that this bulk solution admits two Killing 
spinors and is therefore a BPS state with N = 2 preserved supersymmetry. 

We would like to end this section by stressing a crucial difference between our solution and 
the "conifold" case (see [|l4|] for related discussions). In the former, the S 2 cycles exist only near 
the apex and thus we do not have a 3-form flux at infinity, whereas in the latter, the base of the 
cone supports such cycles and, as a result, the flux at infinity is non-vanishing^]. 

3.2 Proof of bulk supersymmetry 

As usual, in order to investigate the supersymmetry properties of the bosonic solution we have 
found it suffices to consider the supersymmetry transformation of the fermionic fields (the grav- 
itino and the dilatino) and impose that, for a Killing spinor, they vanish identically on the 
chosen background. By using the formulation of |3"2||, one easily gets: 



Hm = V MX + ^iT A ^ F Al ,„ A5 V^T BX 

+h (-^b Ai ...a 3 V* + 9T A ^V^) n +]AlA2A3 X * 
5\ = -ilT A ^n +lAlA2A3 x (3.43) 
where the supersymmetry parameter % is a complex ten-dimensional Weyl spinor: 

rn X = X (3.44) 

and where we have already used the information that on our background the dilaton and the 
Ramond scalar vanish. Using the gamma matrix basis described in appendix |A| (see in particular 
eq.s (A.7,A.£ . A.9D ) and writing the 32-component spinor x as a tensor product: 



X = £ <8> V (3-45) 

of a 4-component spinor e, related to the 3-brane world volume with an 8-component spinor r\ 
related to the transverse manifold Mq, the dilatino transformation rule takes the following form: 

SX = -il 75 e ® T^ k n +lijk v (3.46) 

Hence the first condition a Killing spinor must satisfy comes from the vanishing of the dilatino 
rule and reads as follows: 

u(n)v = r ijk n +]ijkV = o (3.47) 

Using the self-duality relation (|3.37[) we can immediately verify that: 

n(W) = n(W)i(l + ir.) (3.48) 
so that in order to solve eq. Q3.47 ) it suffices to impose the M§ chirality choice: 



ir,r] = -r] (3.49) 



We thank I. Klebanov for clarifying this point. 
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Since Tn = — ij^ ® r,, eq. (|3.49|) implies: 



75 £ = £ 



(3.50) 



Inserting conditions ( |3.49 ) and ( 3.50[) into the gravitino transformation rule, by use of the specific 
structure of Fa 1 ...a$ defined by eq.( 3,33j ), we obtain two new conditions that must be imposed 
on the Killing spinor: 



= V aX ~ Ifila^T 1 X- ^75 la U(H) X * 



(3.51) 
(3.52) 



In the above equations fi is the intrinsic derivative of the harmonic function H defined by eq. 
( p. 34 ). Recalling the form of the spin connection on our background (see eq.fl3.8j)) we easily 
verify that eq. ( 3.51] ) reduces to: 



d a £ ® V ~ 4757a£* ® Tl(H)rf = 



(3.53) 



This equation is satisfied if we choose constant supersymmetry parameters on the 4-dimensional 
world volume: 

d a e = (3.54) 



and imposing the further condition 



on the transverse commuting spinor rf 



U(H) rf = 



(3.55) 



Taking condition ( 3.55 ) into account and recalling once again the form of the spin connection 

= e ® (v + l H~ l dH) i] (3.56) 



(|3.8D, eq. (3.52) becomes: 



where the hatted operators T> and d are the covariant and ordinary exterior derivatives on 
Mq = M 2 x ALE considered by itself without the warp factor H(z, z, y). 

Introducing a rescaling factor H~ l l^(z,z,y) and a further tensor product splitting of the 8 
component spinor 

77 = H- 1/S (z,z,y)9 ® f (3.57) 

into a 2-component spinor 9 relative to the manifold M 2 times a 4-component spinor £ living on 
the ALE manifold, eq. ( 3,56j ) reduces to the condition of covariant constancy on ALE: 



= V£ 



(3.58) 



Collecting all the conditions we have found and using the gamma matrix basis described in eq.s 
(|A.1Q ) of the appendix we can conclude that our background admits as many 4-dimensional 
supersymmetries on the world-volume of the 3 brane as there are simultaneous solutions of the 
following three equations: 



£ = 0-3 9 ® 75 £ 

= vi 



(3.59) 
(3.60) 
(3.61) 



where co^ v are the intrinsic components of the anti-self-dual normalizable 2-forms on ALE. 
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So far we have not used the specific geometric properties of the ALE manifold. It is at 



this stage they come into the game, playing an essential role. Squaring eq.(3.61) we obtain the 
integrability condition 



-\R uv luvt = U (3-62) 

where R uv ^y uv is the curvature 2-form of the ALE manifold. This latter is HyperKahler and as 
such it has a triplet of covariantly constant self-dual 2-forms fP, (x = 1,2,3) whose intrinsic 
components satisfy the quaternionic algebra. This implies that the holonomy of the manifold 
is SU(2) L rather than SO(4) = SU(2)l x SU(2)r and that the curvature 2-form is anti-self- 
dual. This follows from the integrability condition for the covariant constancy of the self-dual 
HyperKahler 2-forms. On the other hand, from the Hirzebruch signature r of the ALE manifold 
it follows that there are exactly r normalizable anti-self-dual forms ujj. From the trivial gamma 
matrix identity: 

75 7™ = -\tuvstl st (3.63) 
follows that the two chirality eigenspaces: 

7 5 £± = H± (3-64) 

are respectively annihilated by the contraction of ^y uv with any self-dual or anti self-dual 2-form. 
Therefore it suffices to choose: 

£ = £+ (3-65) 



and equations ( p.60| ), ( |3.62 ) are simultaneously satisfied. Eq. (3.59) is satisfied by setting: 



a 3 9 = 9 (3.66) 



and once the integrability condition ( |3.62| ) is fulfilled eq. (3.61) can be integrated yielding 



two linear independent solutions £^.' 2 spanning the irreducible representation (0,1/2) of SO(4). 
The other irreducible representation (1/2,0) corresponds to spinors that are not Killing and do 
not generate supersymmetries preserved by the background. In conclusion we have 2 Killing 
spinors generating an M = 2 supersymmetry on the world volume. In other words the bosonic 
background we have constructed corresponds to a BPS state preserving a total of 2 x 4 = 8 
supercharges. 



4 The Eguchi-Hanson case 

As showed in the previous sections, the complete integration of the supergravity field equations 
is reduced to the solution of a single differential equation, namely eq. ( 3.271 ). In this section 



we investigate the properties of such an equation, choosing the simplest instance of an ALE 
manifold, namely the Eguchi-Hanson space , which in the ADE classification corresponds to 
A\ ~ Z2 (see appendix [B] for a more detailed treatment of the Eguchi-Hanson space and for the 
definition of the relevant conventions). While in this section we will go through a somewhat 
detailed mathematical analysis, section 5 contains a summary of the most relevant results and 
their physical implications. It should be said from the beginning that, whereas section 3 is 
entirely based on type IIB supergravity, in the following we are going to make some (we believe 
resonable) assumptions on the nature of the microscopic theory, namely the structure of the 
source terms needed to fix the boundary conditions. 
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The Eguchi-Hanson metric has the form: 



ds 2 EH = g(r)- 2 dr 2 + r 2 (a 2 x + a 2 ) + r 2 g(r) 2 a 2 z 



g(ry = (4.1) 



where r > a, 8 £ [0, tt[ and tp, (ft £ [0, 27r[. The 1-forms cr* = {cr x , a 37 , a 2 } are listed in eq. (B.2) 
of appendix |B|. 

This space has a unique homology 2-cycle S located at r = a and spanned by the coordinates 
9, <p. The anti-self dual form to fulfilling eqs. (3.19) is: 



2 



2 7T V r 

The function A(y) defined in ( 3.1S| ) is explicitly evaluated to be: 



" = — d (-i) (4-2) 



A(y) = 4^ (4.3) 
7r z r 5 

The equation for H in the Eguchi-Hanson case can then be easily obtained from the general 
expression in eq.( |3.27| ). As it is usually the case, we make a spherically symmetric ansatz, 
compatible with the background at hand: for the case we shall be interested in, the coefficients 
of the equation for H depend only on the radial coordinates p = -y/zE and r on M 2 and on the 
Eguchi-Hanson space respectively, and we can then assume the same property for H to hold and 
write it as H(p, r)[|. The equation for H can then be written as: 

(dl + dt) H+±d r (r 3 g(r) 2 d r H) = -J^ \d z 7 | 2 + S(p, r) (4.4) 

where S(p, r) is a source term for the 3-brane charge for which we also make a spherical ansatz. 
Differently from the first term on the right hand side, S(p, r) is not deduced from the dynamics 
in the bulk. In the present analysis its presence should be intended only for the sake of fixing 
boundary conditions near the cycle and it will be determined in the following sections. 



4.1 Solution method by partial Fourier transform 

The method we use to study eq. fl4.4j ) is based on a partial Fourier transform on the coordinates 
of the flat M. 2 space that reduces our problem to that of a massive Laplace equation on the 
residual four-dimensional manifold in the unknown function H(p,r), Fourier transform in M 2 of 
H(p,r) — 1 (with p being the modulo of the momentum in M 2 ). We refer for the details of the 
analysis to appendix [C]. As we are going to show this method allows a parallel treatment of our 
case M 2 x ALEeh with the case M 2 x M 4 which is the covering space for fractional D3 branes 
on orbifolds Q. Such a parallelism is very instructive since the comparison with the flat case 
provides the necessary insight to discuss and interpret the appropriate boundary conditions. 
The main goal of the analysis in the following sub-sections is to determine H(p,r) through 
its near-cycle and far-from-cycle asymptotic limits and an exact analytical description through 
a power series expansion of H in an open neighborhood of the cycle, consistently with the 
boundary conditions. The reason why we choose to study the equation in the Fourier transform 

7 We shall consider the case in which ~/(z) = log (z). In a more general situation \dj(z)\ 2 could also depend on 
the angular coordinate on R 2 . In this case the function H would have an angular dependence in R 2 as well. 
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H(p, r) of H (p, r) — 1 instead of H (p, r) is a matter of convenience in dealing with the boundary 
conditions at infinity: at large distances H (p, r) — ► 1 and therefore H (p, r) — > 0. The relation 
between -ff (^, r) and r) is deduced from the standard definition of Fourier transform given 
in appendix ^ and reads: 

H(p,r) = ±-J™dppJ (pp) (F(p,r)-1) (4.5) 

The equation for H (p, r) has the following form: 

L r [a,p]H(p,r) = J(p,r) (4.6) 
where we have introduced the second order differential operator 

Lr[a,ix]u = ~ ^g(r) 2 ^j-p 2 u (4.7) 

and where we have defined the source function: 

2a 4 

J(/x,r = -^-gj(/x ) + S(p,r) (4.8) 
7r z r° 

j(p) denoting the Fourier transform of |97(z)| 2 and S(p, r) being the transform of the source 
term S(p, r) for the 3-brane charge. 

The notation L r [a, p] is very convenient since it encodes various cases into a single formula. 
For a = 0, L r [0, p] is just the massive radial Laplacian on flat M 4 space. Alternatively for p = 0, 
L r [a, 0] is the massless radial Laplacian on either ALEeh or M depending on the value of a. By 
massive or massless Laplacian we mean the operator (□ — p 2 ) with p 2 > or p 2 = 0. 

It is known from the standard theory of differential equations that the general integral 
H(p, r) of eq. ( |4.6| ) has the form: 

H(p, r) = (3i ui{p, r) + (3 2 u 2 (p, r) + u n _ h (p, r) (4.9) 

where u\2(p, r) are two independent solutions of the homogeneous equation associated with 
&■ 

L r [a,p]u = (4.10) 

while ii n _/i(//, r) denotes a particular solution of the non-homogeneous equation. The constants 
Pi and @2 are fixed by the boundary conditions. 

It is useful to express the unique solution of eq. (|4.6|) fulfilling the required boundary condi- 
tions in terms of a Green function G(a, p\r, r'): 

~ r°° 

H(p,r)= G(a,p\r,r')J (p, r') (r'f dr' (4.11) 

J a 

where the Green function is defined by the equation: 

L r [a, p] G (a, p\r, r') = (4.12) 

The defining properties of G(a, p\r, r') are recalled and thoroughly discussed in appendix |C]. 
The Green function associated with the homogeneous equation ( [4.10 ) is generally expressed in 
terms of a combination of the solutions ui^- The main problem to be dealt with in the next 
sections is to determine, on one hand, the source term S(p, r), and, on the other hand, the 
form of the Green function such that it implements the required physical boundary conditions 
of the problem. These are conditions on the behavior of H near the cycle r — > a and at infinity, 
r — > oo. 
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4.2 Asymptotics of the the two independent solutions ui,u 2 for the homoge- 
neous equation 

So far our discussion of the differential equation has been fully general and applies to both R 4 and 
to the Eguchi-Hanson case. We turn next to a more detailed consideration of the homogeneous 
equation and of its basis of solutions in various sub-cases. Indeed the parameter space of our 
differential operator, spanned by a and fi is partitioned into four sub-cases according to the 
following scheme: 





a = 


a / 


fi = 


massless R 4 case 


massless EH case 




massive R 


massive EH case 



In three of the above four sectors the homogeneous equation can be solved in terms of elementary 
transcendental or special functions so that all asymptotic behaviors are easily calculated. The 
really hard corner is the lowest right one, namely the massive Eguchi-Hanson case. As we show 



in the next section, for /x ^ 0,a / the homogeneous equation (4.6) is a particular instance of 
the confluent Heun differential equation and its solutions can be expressed in terms of a power 
series. The same pattern of asymptotic behavior is reproduced also in this case and provides 
the needed guideline to fix the boundary conditions. 

Let us write down the basis of solutions in the three cases where they are easily accessible 
by elementary methods and indicate our nomenclature for the massive EH case that we address 
in the following section: 



a = Q;[A = massless R 4 ui(0,0, r 

u 2 (0,0,r 

a = 0;/j,^0; massive R 4 ui(0, fi,r 

u 2 (0,^,r 

a^0;/x = massless EH u\(a,0,r 

u 2 (a,0,r 

a^0;/i/0 massive EH ui(a,[i,r 



1 



•2^ (/) (/V 

H( R \n,r) 



(4.13) 



where K u (z) and I v {z) denote Bessel function of index u and type K or I, respectively. The 
reader can easily check that the above solutions have been normalized in such a way as to 
fulfill the normalization ( C.18| ) of the Wronskian. The notation iT^O anc l _ff( R ) is reserved for 
the confluent Heun functions, to be discussed in the next section, that have the corresponding 
behavior at r — > a as their analogues in the massless case. To understand the rationale of the 
previous statement it is instructive to consider the series expansion of 1^2(0, /J,, r) at r = 0. We 
find: 



ui(0,/x,r) ~ r ^ 



4 + ^ 



-l + 2 7 + 21og 



v 2 



+ 0{r 2 ) 



2 2 



(4.14) 



(4.15) 
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where 7 = 0.577216 is the Euler gamma number. As one sees the leading behavior of the massive 
solutions for r — > is just identical with the behavior of the massless solutions, u\ is singular and 
has a second order pole 1/r 2 while 112 is a regular power series and tends to 1. Indeed choosing 
the overall normalization consistent with eq. ( C.18| ) we could also arrange for the coefficients of 
the pole and of the constant term in u\ and K2 respectively, to be equal to one. In the case of 
the regular solution U2(0,fi,r) the generic coefficient of the power series expansion in eq.( 4.15 ) 
can be easily calculated and we have: 

«a(0, fi, r) = 2 4 n(( 1+ ")i)2 (^ r ) 2n ( 4 - 16 ) 

As we see in the next section a similar regular series, although much more complicated can be 
written down for 112(0, fi,r) = h re9 (a, fi,r). Let us next consider the behavior at r = a of the 
massless solution U2(a, 0, r). We have 



m(a,0,r) ~ r ^ a log \+0{r-a) (4.17) 

2 cr \ a J 

As we see the double pole singularity 1/r 2 appearing in the M case has been replaced by a 
logarithmic singularity log(r — a) appearing in the Eguchi-Hanson case. What remains identical 
in the two cases is the relation between the massive and massless solutions. Also for // / we 
find in the next section that: 

log(r 
~2^2 

u 2 {a,fi,r) ~ r ^ l + 0(r-a) (4.18) 



ui(a,fi,r) ~ r _ a — hreg 

1 n & 



Using this information on the asymptotic behavior of the two fundamental solutions we can 
complete our program of constructing the Green function by imposing the boundary conditions 
of physical relevance. 

To this effect we need the preliminary analysis of the near-cycle limit in the Eguchi-Hanson 
space. Let us go back to the Eguchi-Hanson metric in the form (4.1) and let us perform the 
following change of variable: 

Vr 4 - a 4 fl . 

r^v = 4.19 

2r v ; 

The form of the metric near the cycle (r — > a) can now be calculated by expanding the metric 

( |4. ip in power series of v at v = 0. We obtain: 

n 2 

ds 2 EH dv 2 + v 2 dip 2 + — (d9 2 + d(j) 2 sin 2 9) + O (v 4 ) (4.20) 

showing that near the homology cycle the Eguchi-Hanson metric approximates that of a manifold 

m 2 x s 2 . 

In view of this result we consider the transcription of the homogeneous radial equation to 
the new variable v, defined by eq.( [4.19D . We obtain, for small v: 

^H+-^-H-li 2 H = (4.21) 

dv z v dv 

admitting the following two linearly independent solutions: 

ui(v,fi) = K (v(i) ; U2(v,(i) = I (v fi) (4.22) 
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whose behavior for v — > is: 



lo g( u ) 

ui{v,n) ~ v _+ hreg 



2 2 

«2(«,/i) ^->o l + ^ + 0(v 3 ) (4.23) 



If in eq.s (4.23) we reinsert the original variable r via the relation (4.19) we obtain the announced 



proof that, near the cycle, the leading behavior of the massive solutions ui(a,fi,r) is the same 
as the behavior of the corresponding massless solutions. By such a substitution we get: 



ui(a,/i,r) ~ r _> m( ^ ( g a \ v) = -2^1og(r - a) + reg 

__ J(r-a) l 4 - Z4 J 

u 2 (a,n,r) ~ r _^ a «2(-^ — ,M) = l + C>(r-a) 

4.3 The 3-brane charge, the physical boundary conditions and the final form 
of the Green function 

Equipped with the results of the previous section, let us now consider the case of a 3-brane with 
vanishing 3-form placed either at the origin of M 6 in the orbifold case where the transverse space 
is Mq = R 2 x M 4 /r or at the homology cycle r = a in the case where the transverse space is 
Mq = M? x ALEeh- Naming Q the charge of such a brane, the expected behavior of the H 
function near the brane is, in the two cases, the following one: 

H ~ — = =- forM 6 ; H ~ -Qfi?-— for M 2 x ALE (4.25) 

(xl + xj + ... + xj) 2 {p 2 + v 2 ) 

The first of eq.s (|4.25|) is obvious. The second is due to the discussion of the previous section. 
Near the cycle the Laplacian on Mg becomes: 

□m 6 ~ n R 4 + D 5 2 (4.26) 

and since S 2 is a compact, positive curvature manifold there are no zero modes of dga except 
the constant. Therefore the non trivial part of H behaves as a harmonic function on M 4 , namely 
the second of eq.s ( [4.25 ). The scale factor l/(a 2 ) appearing there is understood as follows. If 



Q is the total charge perceived at infinity, the density of charge on the homology 2-sphere of 
radius | is: 

CT7T 

Hence what appears as charge in the near cycle M 4 plane is q rather than Q. Finally the tt/2 
factor that is needed to match the second with the first of equations ( |4,25| ) is just a matter of 
convenient normalization. Let us now perform the Fourier-Bessel transform of eq.s ( [1.25 ). We 
obtain 

1 f° t ( \ 9. a Q V xr ( \ Q tn \ 

Wo Jo(/ip)p FT7F 9 = 4^7 Kl(/xr) = 4^ ni(0 ^' r) 

hL Mfl(>)p 2a 2 (v 2 + p 2 ) dp = ^h Koi/iv) = ^ 1IlM (428) 

The important conclusion implied by the above analysis is that, in terms of the fundamental 
basis of solutions u\,U2 we have obtained the physically appropriate boundary conditions for 
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the function H(p,,r). In both the orbifold or smooth ALE case, in the limits r — > and r — ► o 
respectively, we have: 

F(//,r) ~ £uf (/x,r)+reg (4.29) 

where Q is the 3-brane charge and u^ w (/i, r) denotes the divergent part of the solution u\(a, fi, r) 
defined in eq.s( [4.13" ). This condition, together with the boundary condition at infinity, fixes the 



coefficients f3\, 02 in eq. (|4.9|) to be: 



Pi = §; ; & = (4.30) 



The appropriate expression for the Green function in order for it to implement the physical 
boundary conditions defined above is derived in appendix ^ and is the following: 

G(r,r )-\_i Ul(T)u2{jl) forr >r' ^.31) 
Acting with such a Green function on a source: 

Sfcr) = 4*^ (4.32) 



we obtain a solution with the prescribed behavior as in eq. ( [4.29| ). From the physical viewpoint 
the source ( [4.32| ) comes from the world-volume action of the 3-brane. 

4.4 Reduction to the confluent Heun equation 

As we have seen in previous sections a convenient approach to the solution of equation fl4.4|) 
relies on the partial Fourier transform of eq.(C.l) leading to the new eq.s ( |4.6| ) , ( [4.7D for the 



function H(/j,r). So far we have analyzed such an equation from a general viewpoint and we 
have discussed how to impose the correct boundary conditions from the asymptotic behavior 
of the two solutions u\ and v>2 of the associated homogeneous equation. In this section we 
eventually come to a more in depth analysis of the equation in a generic point of its moduli 
space, i.e. at ^ ^ 0, a ^ 0. 

Parameterizing the radial direction in the Eguchi-Hanson space through the new variable 
w = (a 2 — r 2 )/(2a 2 ), equation ( [4,6D ,( |i~7|) can be rewritten in the following form: 

o d 2 H , . dH ~ h{n) w 



dw 2 dw ' (I — w)(l — 2w) 3 

H{H,w) w ^°° (4.33) 



where 



(4.34) 



We have omitted in the source term the part proportional to S(fJ,, r) since, as explained in 
the previous sections, its effect amounts in a determination of the relative coefficients of the 



two independent solutions u\ 2 of the homogeneous equation, as in eq. ((4.30|). We postpone to 
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appendix a more detailed analysis of eq. (|4.33[) , as a non- homogeneous confluent Heun equation, 
and of its solution. In the present section we shall determine the analytical expression of H in a 
finite neighborhood of the cycle, represented by the regular singularity w = and implement on 
it the the required boundary conditions: the behavior of the integral at infinity, defined in eq. 
( [4.33 ), and its behavior at w — > 0, determined by the presence of source terms as we explained 



in section [4^. Finally, using these exact results we shall compute the anti-Fourier transform H 
of H in its asymptotic near-cycle limit, confirming the predictions made in the previous section. 

Let us now discuss the solution around w = and the asymptotic solution for uu — > — oo. 
The expansion of the solution around w = corresponds to considering its behavior near the 
exceptional cycle: r ~ a. Since w = is a regular singularity there are standard ways to 



determine the power series expansion of the the solutions for both the homogeneous [37] and 
non-homogeneous equations around such a point. 



First let us focus on the homogeneous equation associated with (4.33) and let us denote it 
symbolically by L w H = 0. This equation has two independent solutions, one being regular at 
w = will be denoted by H^(p, w), the other having a logarithmic behavior near the same 
point will be denoted by H^ J \^, uu). From standard analysis we may determine the power series 
expansions of the above integrals around w = 0: 

oo oo 
n=0 n=0 

H^(p, uu) = lim 4- I w ° Y c n(^, °)w n ] (4.35) 
da- \ ^ J 

the coefficients c n (/i, a) will be derived in appendix ^] and have the form: 



F n {ji,&) = detF n 0i,a) (4.36) 
where F n (/i, a) is an n x n matrix, see appendix |C], whose entries are: 

PMas <r)]u = k-a-n + i; [F„(/i, <r)] fi>i = -{k + a + n - j) 

<r)] i+u = (<J + n-i) 2 ; i,j = l,...,n (4.37) 

Both the series for H^(fi, w) and H^ R \ii, w) converge for \w\ < 1, namely for r 2 — a 2 < 2a 2 . 
A generic solution of the homogeneous equation will be expressed as a linear combination of 
H^\fi, w) and H^ R \fi, w), the relative coefficients will be determined in quite the same way 



as the constants @\ , defined by eq. ( p^Gj ) were determined in section [O 



Let us now consider a particular series solution of the full non-homogeneous equation 

oo 

H {n ~ h) (/i, w) = J2 <~ h ™ n ( 4 - 38 ) 

n=0 

the coefficients c^~ h are derived in appendix [C] and have the form 



nj=i/o(7) n-=ii 2 

F n (/i) = detF n (/i) (4.39) 
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where F„(/i) is an n x n matrix similar in form to F n (/i, a) except for the last column which 
depends on the source term expansion: 

PM/i, <t)]u = k-a-n + i; [F n ( M , a)] ij>{ = -(k + a + n - j) 
Wn(p,<r)] i+ u = (cr + n-i) 2 ; 

Wn(p, cr)] in = -hit 1 ) s n-i ! z = 1, . . . , n ; j = 1, . . . , n — 1 (4.40) 

where the coefficients Sj are defined by the power series expansion of the source term on the 
right hand side of eq. fl4.33| ): 



Kn)w , M f 



+ 1 



(1-^(1-2^)3 



^(i + l)(i + 2)2- 1 (4.41) 



i=0 



it can be checked that both the series expansions for the source term and for H^ n ~ h \p, w) 
converge for < 1/2, i.e. for r 2 — a 2 < a 2 . 

Below we write the expansions of H^ifi, w), H^in, w) and H^ n ~ h \p, w) up to order u; 3 : 

fa w ) = c Rfa _ kw + ^ w 2 + 0(u> 3 )^) 

H^faw) = c^fa (log (w) + (-k log (w) + l + 2k) w 
+ ^ (A: 2 log (to) + 2 - 3A: 2 ) w 2 + . . ^ 

H^'^faw) = hill) (w + ^(S-k)w 2 + 0(w 3 )^j (4.42) 

In some particular limits the above series can be summed. 

First of all let us consider the solution at p 3> a to which only the Fourier transform H at 
p = (k = 0) contributes. In this case, the series in eqs. ( [4.42 ) are easily summed and we get: 

H^fa w) = 1 (4.43) 

H^'in, w) = log (w) + — = log (—w) — log (1 — w) + const. 



n=l 

(r 2 - a 2 \ 

= log -= 5- + const. (4.44) 

\r z + a z / 

H^^,w) = - ^ = + const. (4.45) 

VP ' ; 2(1 - 2w) 2 2r 2 v ' 

Eq.s (|44l) and flOg) are a very nice confirmation of what we had already established in previous 
sections considering the limiting form of the equation either for /i — > or for a — > 0. In particular 
eq.s ( [4.43 ) and ( [4.44 ) justify the absolute normalization of the solutions iti^ chosen in eq.( 4.13[ ). 

Another accessible limit is that corresponding to /x 2 3> 1/a 2 with \w\ *C 1 (i.e. r 2 — a 2 <C a 2 ) 
and keeping < 1. This is the asymptotic near-cycle limit in which we approach the brane 
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also along M 2 . We can sum the series up to 0(1 /k 2 ) for the solutions of the homogeneous 
equation and up to 0(1/ A; 3 ) for the particular solution of the non-homogeneous one: 



w 



W 



Hi R \kw) + \Hf\k W )+0{^) 

H { I \kw) + ^H[ I \kw)+O(^) 
1 - T (n-/i)'- x 1 



r (n-h) 



k 



{kw)+0 {-) 



(4.46) 



The functions H[(kw) in the above expressions are obtained by summing the terms propor- 
tional to k n ~ l w n for the solutions of the homogeneous equation and the terms proportional to 
fcn-l-i w n - m foe series for H^ n ~ h \p, w). This requires the computation as a function of n of 
the coefficient of the k n ~ 1 or k n ~ l ~ 1 order term in the Hessenberg determinants F n (/j,,0) and 
F n (n). Such a computation can be done but, for simplicity, we focus only on the cases 1 = 0, 1). 
It can be easily checkeof] that for I = we have: 



H, 



o 

(n-h) 




(0 



£71 

HfHO = ^^^ = ^10(2^) 

n=0 

-h(n)I (2y/£)+h(jJ.) 



(4.47) 



where we have set 



\kw\ 



2 i 2 



a 2 )/8 < 1 



(4.48) 



Once again, eq.s (4.47) are a nice confirmation of what we had predicted in the previous section 
by the analysis of the differential equation in the near cycle limit. ^From our general power 
series solutions we retrieve the result that the regular and irregular solutions of the homogeneous 
equation in eqs. ( [4.47 ) have the same form as the massive Green functions on M 2 . For a detailed 
comparison with our previous results it suffices to observe that: 



jl V = jJL- 



2r 



\fr 2 



V2 



(4.49) 



Eq.s fl4.47| ) allow also to identify the correct boundary values of the constants appearing in the 
homogeneous solution. As announced we must have: 



Q 

Ana 2 







(4.50) 



Indeed eq.s ( |4.47| ) show that near the cycle the particular integral H( n ~ h ^ goes to zero as foreseen 
in our general discussion. 

At first order in \jk a particular integral of the non- homogeneous equation, is simply 
w) ~ h(fi)/k and therefore the behavior of its Fourier transform H^ n h \p, r) in 
the near-cycle limit is : 



H {n - h \p, 



r) ~ 



2tt 2 a 4 



Po 



(4.51) 



8 For simplicity we assume that the coefficients Cq and Cq do not depend on , 
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where according to the discussion in section 3 we have assumed the holomorphic function 7 to 
be a logarithm (see also the discussion in appendix C), and po is an integration constant acting 
as a regulator, not to be confused with the enhangon radius. The next order terms (I = 1) in 
eq.s (4.46) are written in appendix^. 



As far as the asymptotic behavior of the solution at infinity r 3> a is concerned we observe 
that the first term in the r.h.s. of the equation for H(p, r), eq. (|4.4| ), is just a resolved delta 
function for radially dependent test functions <3?(r) in the ALEeh space: 



lim 



/ 


' 4 


a 4 " 


1 ALE 




^.8 



Q(r)dn ALE = $(0) (4.52) 



Therefore for r 3> a the differential equation (|4.4j ) becomes homogeneous, namely its r.h.s. 
vanishes. The corresponding equation in the Fourier transform H(p, r), for w — > —00 has two 
independent solutions of which only one is consistent with the boundary condition at infinity 
( [4.33 ), namely: 



H(p, r) = c — Ki(fir) ~ c X ui(a, fj,, r) (4.53) 
r 

where we have used the radial variable r instead of w ~ — r 2 /(2a 2 ). The coefficient c has already 
been fixed in terms of the 3-brane charge and the absence of the near cycle regular solution, 
that is instead divergent at infinity is confirmed. As we expect in this limit, the behavior of the 
warp factor H is just the same as in the orbifold case (7j, schematically: 

(p, r > a) : H(p, r) ~ ® (1 + log terms) (4.54) 

On the other hand, in the region near the cycle, combining the asymptotic behaviors of the 
homogenoeus and inhomogeneous solutions, we find 

P<a ): H(p,r)~ K log (P-X (4.55) 

r~aj VP ' ; 2a 2 p 2 + v 2 2tt 2 a 4 \PoJ 

where if is a (undetermined at the level of the present analysis) constant. The corrections to the 
above behavior can be systematically deduced from the power series expansion of the solution 
described in the present section. 

Looking at eq.(4.55), it is clear that there could be a value of r, p for which H = 0, this being 



an indication for the presence of a naked singularity of the repulson type ||21[| . This singularity 
should be removed, somehow. In most non-confomal M = 2 versions of the gauge/gravity corre- 
spondence this singularity has been shown to be excised by the so-called enhangon mechanism 
p3| l. This is the case, for instance, of fractional branes on orbifolds, 0, [jO]]. The value of the 
enhangon corresponds to p = p e , the scale where the scalar field 7 vanishes. This in general 
turns out to be the scale where the dual gauge theory becomes strongly coupled and new light 
degrees of freedom are expected to become relevant, both at the gauge theory (where instanton 
effects become important) and at the supergravity level (where tensionless strings occur). For 
all this analysis to work, it is important that p e is bigger than the scale at which the repulson 
occurs. In fact, when this is the case, the region where supergravity is reliable, namely p > p e , 
is free of any singularity. In order to see if this happens also in our case and if the cut-off p e has 
in fact the expected meaning, one should have a full control on the world-volume action of the 
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source, with all numerical coefficient, normalization, etc... This we cannot do at the level of the 
present pure supergravity analysis. 

Our solution, however, differs from that of fractional branes on singular space because of 
a better behaviour of the warp factor H on the plane M? where it is non-singular, while the 
solution of the field 7, which is responsible for the enhangon mechanism, has essentially the 
same structure. It is then natural to believe that the nice properties described above should not 
drastically change. Hence we consider our solution to be reliable, well-defined and singularity 
free for p > p e (where the leading order near cycle behaviour is well described by the first term 
in eq. (|4.55|) itself) while we postpone a more precise analysis of this issue to future work. 



5 Discussion and Conclusions 

In this paper we have constructed a new class of solutions of type IIB supergravity describing 
3-branes whose transverse geometry is of the type Mq = M 2 x ALE for an arbitrary smooth ALE 
manifold. We have shown that these solutions preserve 8 supercharges and thus correspond to 
M = 2 supersymmetry in four dimensions. They are characterized by constant dilaton and axion 
field, a 3-brane metric 

ds 2 = {-ri^dxf* dx v ) + mdzdz + H^ds 2 ALE , (5.1) 

a self dual five-form 

F^f = d (H- 1 fi Rl , 3 ) + *d (H^ 1 fi Rll3 ) , (5.2) 

and a complex three-form, valued only in the transverse space, whose form is uniquely fixed by 
supersymmetry: 

n + = 2dr/i(z) A uj 1 (5.3) 

(7/(21) being analytic functions on M 2 and loj harmonic square-integrable anti-self-dual forms on 
ALE). The whole set of type IIB equations of motion reduces to solving the following PDE for 
the warp factor H 

(□ R2 + Dale) H = -d zlI drfj A IJ (y) (5.4) 

This equation is quite more involved compared with the analogous equations for the warp factor 
arising from N = 1 supersymmetry which are of the ODE type due to the conical nature of 



the transverse space. We have then proceeded in studying eq. (5.4) in the simplest instance 
of an ALE manifold, the Eguchi-Hanson space, whose properties were reviewed in section [| 
and appendix |B|. In the case in which the (here unique) analytic function 7(2) possesses radial 
symmetry, the equation for the warp factor took the form 



±d pP d p H+±d r (r 3 g(r) 2 d r H)=-^\d 1 (z)\ 2 , (5.5) 

and, after transforming to the variable p conjugate to p, reduces to a particular type of confluent 
Heun equation, whose solution by power series can be determined with arbitrary accuracy. 
Particularly interesting is the behavior of H(p, r) near the cycle r = a. It is still necessary to 



invoke an enhangon-like mechanism [23] to cut-off the region p < p e where supergravity is not 
reliable. However, contrary to the orbifold case, the warp factor H(p, r) is regular as r — ► a 
for p > p e , due to the fact that the contribution of the three- form to the source term has been 
regularized. In other words, wherever supergravity is applicable, smoothing out the transverse 
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geometry does regularize the solution for the warp factor. The above result can be understood 
by studying the limit of eq. (|5.5| ) in the presence of an auxiliary source term 

source oc -Q 6 ^ ~ "V (z) (5.6) 

on the r.h.s. introduced to fix the boundary conditions near the cycle. In the region r — > oo 
the Laplacian reduces to the flat Laplacian on R 6 and Q is the total charge observed at infinity 
(provided the complex 3-from is vanishing). In the opposite limit, r — > a, the Laplacian reduces 
to the flat Laplacian on M 4 x S 2 (the S 2 dependence of H is trivial) and the sources are clearly 
regular as long as p / 0, thus so is the solution. By computing the density of charge q on 
the two-cycle one finds that it is consistent with the interpretation that the three-brane charge 
is uniformly smeared over the cycle (q = Q/na 2 , being ira 2 the area of the cycle), thus each 
element cE of the surface area of the cycle contributes q cG to the charge at infinity. The simple 
form of the Laplacian near the two-cycle, and the fact that we have control over the boundary 
conditions, allows us to write down the asymptotic form of the metric in that region. Neglecting 
the contribution of the three-form compared with the source term for the five-form [] we find 



d p pd p + -d v vd v H cx -q5(v 2 )5(p z ) (5.7) 
p v J 

which is nothing but the Laplace equation for a point charge in four dimensions written in "bi- 
polar" coordinates. Introducing the variable R = \J (p 2 + v 2 ), denoting the radial coordinate of 
R 4 in units of q we find that H ~ q/R 2 , yielding 

ds 2 ~ 4= (-V^dx^dx") + ^-dR 2 + ^Rds 2 s3 + ^-^yds 2 s2 (5.8) 
y/q R R 4 

where ds| 2 is the metric of the two cycle of Eguchi-Hanson and ds| 3 is the metric of the three- 
sphere at fixed R in M 4 . An interesting point is that the above result holds even in the absence 
of the three-form flux and shows that conformal invariance of the dual theory is always broken 
because the metric in R and x^ is no longer of anti de Sitter type. Obviously we cannot take the 
limit R — > 0, because this would imply going to large curvatures, and in particular crossing the 
enhangon radius. However, the result suggests that a possible interpretation for the parameter 
a in the Eguchi-Hanson metric is that of a Fayet Iliopoulos term, breaking conformal invariance 
in the infrared, in accord with previous work on the subject. 

The situation can be summarized as follows. In the absence of flux the exact supersymmetric 
3-brane solution that we have found interpolates between a standard ten dimensional 3-brane 
solution at the singularity of the metric cone on S 5 , i.e. the standard R 6 manifold, and a 3- 



brane solution of an effective 8-dimensional supergravity. As explained in the literature [38] 
(for a review see |^9|) we can always consider sphere reductions of all supergravity theories and 
in particular an S 2 -reduction of type IIB supergravity. This yields an effective 8-dimensional 
supergravity that has 3-brane solutions. In this case however, there is a coupling to an effective 
dilaton that emerges as the conformal factor of the metric in the dimensional reduction. Hence 
the 3-brane solution of the 8-dimensional supergravity is no longer conformal and we can follow 



the prescription of Townsend et al [40] by making the transition to a dual frame where the metric 



9 The argument below becomes exact in the limit where the three-form charge vanishes and it is applicable as 
long as the p e cut-off is not too small. 
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factorizes into the product of a 3-sphere metric times a domain wall solution of an effective five- 
dimensional supergravity theory. 

Although the main thrust of this paper has been on the supergravity solution, there are 
various properties of the dual gauge theory that can be described to complete the analysis. 

We have seen that supergravity alone is not sufficient to determine 7(2:) or the boundary 
conditions for the warp factor, while this is easy in the case of the orbifold (by using for example 
boundary state techniques) or on the conifold (in which one knows a priori the structure of the 
source term, despite one does not know the world-volume action). This arbitrariness implies 
that our solution describes various deformation or various vacua of M = 2 theories. 

Let us start by taking the orbifold limit and 7(2) = const. What we obtain is clearly (N) 
regular D3-branes at the orbifold point. The constituent branes could move away from the 
orbifold point in all six transverse directions, thus the dual gauge theory has the relevant hyper- 
multiplets and it is the well known conformal version of an N = 2 Yang Mills with product 
gauge group 11/= 1 U(di N) (k being the number of irrepses of the discrete kleinian group V C 
SU(2) that characterizes the orbifold and di their dimension) and hyper-multiplets in the bi- 
fundamental representation. The blow up constitutes no problem in this case, and it is well 
understood in terms of Fayet-Iliopoulos terms that break conformal invariance but preserve 
super symmetry. In the simplest case, namely for r = Z2, which is the one we focused on, there 
are three such terms. 

On the other hand, in the orbifold limit (a —* 0), it is known that we may switch on a 
three form flux described by j(z) 7^ const, which preserves the N = 2 supersymmetry. In the 
unbroken phase symmetry requires that we can only have j(z) = {£ > 0) or 7(2) = log (z) 
(£ = 0). In the case of a small flux, the corresponding operator on the boundary CFT side can 
be seen as a perturbation with conformal weight A = £ + 4. This perturbation is irrelevant 
for £ > 0, not affecting therefore the conformal IR fixed point, and it is marginal for £ = 0, 
signaling the presence of a small number of fractional branes together with the regular ones at 
the orbifold point. 

The original and perhaps surprising feature of our solution consists in the fact that it may be 
thought of as the simultaneous effect on the CFT describing regular D3~branes at the orbifold 
point of the two kind of deformations discussed above: the relevant one related to the blow-up 
moduli and the marginal/irrelevant one related to the non-trivial flux. Indeed, as previously 
pointed out, the solution j(z) ~ log (z) occurs in the orbifold case for the fractional D3-brane 
solution where the need for a logarithm is dictated by the source term coming from the world 
volume action |/J . Yet we know that in that case the orbifold cannot be blown up [^4|] and this 
can also be understood from the fact that the gauge theory dual to the fractional brane lacks the 
fields necessary to construct the appropriate Fayet-Iliopoulos terms. One might have deemed 
that this fact would rule out solutions of this kind in type IIB supergravity: on the contrary, in 
this paper we have shown that such solutions are indeed possible, although they might not be 
interpreted as the resolution of fractional branes on orbifolds. 

For small fluxes the most obvious interpretation is that we are perturbing the (already 
resolved) theory on the regular brane as discussed above, but there is also the possibility that 
some backgrounds with logarithmic behavior correspond to blowing up a fractional brane not 
at the orbifold point (background Bjys field flux through the cycle equal to 1/2 in string units) 
but at the singular point (with zero flux), where the no-go theorems [^] do not apply. This 
raises the interesting issue of whether our solution can, in some sense, be seen as interpolating 
between the SCFT + FI theory with product gauge group and the pure YM theory typical of a 
fractional brane. 
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This is one instance in which supergravity alone is insufficient in determining which is the 
dual gauge theory and we need string theory inputs as discussed above. This is not too surprising 
in the case of N = 2 theories since, strictly speaking, supergravity alone is never sufficient to 
determine a dual theory a la Maldacena, (for instance we cannot take the near horizon limit 
without invoking the enhangon mechanism). 

In the future, we plan to return to these issues by making a careful study of the boundary 
action and trying to test if the above supergravity picture can be confirmed by a microscopic 
analysis. It would also be interesting to make connection with very recent work on D-instantons 
M] and wrapped fivebranes [42] and [43]. 
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A Notations and Conventions 

The supersymmetric formulation of type IIB supergravity we rely on is that of Castellani and 



Pesando [[3^] that uses the rheonomy approach [33]. Hence, as it is customary in all the rheonomy 
constructions, the adopted signature of the space-time metric dr 2 used to set the conventions 
for the gamma matrices is the mostly minus signature: 



rjAB = diag < 




(A.I) 



However, in presenting the solution for the metric we used ds 2 = —dr 2 to conform with more 
recent usage. 

The index conventions are the following ones: 

A,B,... = 0, • • • , 9 Flat indices in D = 10 

a,b, . . . = 0, • • • , 3 Flat indices on M 1 ' 3 

a, f3, . . . = 4, 5 Flat indices on IR 2 

t,s, . . . = 6, 7, 8, 9 Flat indices on ALE 

= 4, . . . , 9 Flat indices on M? x ALE) 

M, N, . . . = Curved indices in D = 10 



Curved indices on R 1,3 



z,z = Curved indices on 
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r, a, . . . = Curved indices on ALE 

I, J,... = 1,2, ... ,k indices labeling ALE harmonic self dual 2-forms 

(A.2) 

The complete vielbein of 10-dimensional space-time are denoted V A , those of the M 2 x ALE 
manifold without warping factor are denoted e l . The convention for constructing the dual of an 
n-form uj in 10-dimensions is the following: 

w = u Al ... An V Al A ... A V An e> *cu = ^ e Al ... Al0 _ nBl ... Bn u; B ^ B " V Al A ... A V A ™~" 

(A.3) 

Note that we also use n-form components with strength one: uj = uJ Al ... An V Al A ... A V An and 
not with strength n\ as it would be the case if we were to write u = h^A\...A n V A ... A V An . 
When it is more appropriate to use curved rather than flat indices then the convention for Hodge 
duality is summarized by the formula: 



* {dx M - A . . . dx M -) = Y 1Q n v g MlJVl • • • 9 MnNn e Lx ... Lw _ nNl ... Nn dx* A... dx L ^ (A.4) 

Applying the rule ( |A.4| ) to the case of the specific 10-dimensional metric of our ansatz ( |3.2[ ) we 
obtain the following useful formulae for Hodge duals: 



* (dz A S7 K i,3) 

★ (dz A £1 A le) 

-k(dy T An R i,s) 
★ (dz A lo 1 ) 

^ALE 



-\H 2 dz A £1 A le 
\H~ 2 dz A R i,3 

~^H 2 (idy T ) Adz Adz = —H 2 (My T ) A dx 4 A dx 5 

— ii7 1 dz A uj 1 A K i,3 
dx° A ... dx 3 

det (g ALE ) dy 1 A... dy 4 



(A.5) 



where * denotes the Hodge dual in the ALE space. The Laplacians on R 2 and the ALE space 
are: 



Dr 2 

a ALE (•) 



o 



dx 5 



d y 

d(ALE)* d{ALE) (■ 

1 8 T 



y/det (g ALE 

For the gamma matrices our conventions are as follows: 



det (g ALE ) (g ALE ) T(7 d a 



(A.6) 



{T A ,T B } 



-AB 



pa 

r 



7 a <8) 1 
75 &> t % 



(A.7) 
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where: 

{r,l b } = 2rj ab ; 75 = -i 7 V7 2 7 3 (Ag) 

are the gamma matrices in Lorentzian four space and on the six dimensional manifold Mq 
respectively. The matrices 75 that anticommutes with all 7° and r, that anticommutes with all 
r* are respectively hermitian and antisymmetric and we have: 

7 5 2 = 1 ; r 2 = -l (A.9) 

Since the signature of Mg is completely negative the t % matrices that realize the associated 
Clifford algebra can be taken to be all real and antisymmetric. They can be further decomposed 
with respect to the submanifolds R 2 and ALE as it follows: 

r l = { . _ (A.10) 

[ t u = 1a 6 <g> 7u 

where a a = a 1,2 and a 3 are the standard Pauli matrices while: 

{7ui7«} = *u» ( A - n ) 

are 4x4 hermitian matrices forming an Euclidean realization of the four-dimensional Clifford 
algebra. 



B The Eguchi-Hanson space 

In the present appendix we state our conventions for the description of the Eguchi-Hanson space. 
The metric is defined in eq. ( |4.1| ) of section ||, where the 1-forms a 1 = {a x a y a z } are an explicit 
realization of the Maurer Cartan forms of the SU(2)-Lie algebra: 

da* = e ijk a j Aa k (B.l) 

and their expression is: 

a x = — — (cos tp d9 + sin 9 sin tp d<j>) 
(j y = — (sin ipdO — sin 8 cos ip d(f>) 

a z = -- (dip + cos 9 dc/)) (B.2) 

Let us now consider the properties of the homology 2-cycle which, in this case, is unique. In the 
coordinates (ip,r,9,<p) the cycle is spanned by (9, <p) setting r = a. At r = a the Eguchi-Hanson 
metric has a Bolt singularity which is a coordinate singularity given that ip has period 2ir. The 
limit r — > a is therefore regular for the solution, as can be verified by changing coordinates from 
r to u = r 2 g(r) 2 and taking u — » 0. Restricting the metric to 9, <p with r = a yields: 




(B.3) 



that describes the exceptional cycle S as a sphere of radius a/2. 
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Next let us write the explicit form of the harmonic anti self-dual 2-form u fulfilling properties 
( p. 19 ). Since a z is not globally defined on £ the form: 



a 2 , /cr z \ a? , , , a 2 „ , , , a 2 



lo = — d( — ) = ^drAdip-\ ? cos Odr Add + - — ^ sin 9 d9 A d<p (B.4) 

2 7r V H / 2irr 6 2irr 6 Airr z 

is closed but not exact and constitutes the only normalizable representative of the second co- 
homology group of the space. It can be checked that: J s to = 1, J ALE to A lo = —1/2 = — C~ l 
and that -ku = —lo. To check the last property we fix our conventions for the Levi-Civita tensor 
etpr6(f> = 1 and give the following useful formulae: 

2 cos 9 2 

* (dO A d(j)) = - dr A dtp dr A dip 

r sin 9 r sin 9 

*(drAdiP) = r — (l- (-Y cos9 2 ) d9 A d<t> - g(r) 2 d9 A dtp 

2sm9 \ \rJ J 2sm.9 

*(drAd<p) = —^—g(r) 2 d9Ad'tp+ r -^-g{r) 2 d9Ad4> (B.5) 
2 sin 9 2 sin 9 

where we have used the expression for the determinant of the Eguchi-Hanson metric: 

det(g ALE ) = (r 3 sin0/8) 2 (B.6) 



Finally the expression ( |4.3| ) for the function A(y) is easily derived from the following wedge 
product: 



r 3 



to A lo = -A(y) VLale ; &ALE = — sin 9 dip A dr A d9 A d(p (B.7) 

8 

C Analysis of the equation for H in the Eguchi-Hanson space 

The present appendix is to put on a more rigorous mathematical ground the results on the study 
of the warp factor H discussed in section || and to fix the relevant conventions. The general 
definition of the Fourier transform H(p,y) (p here is the momentum vector along M 2 ) is 

~,_ x i I dz Adz , , . 

H(P,V) = 27 ~(2^f~ exPl^ + ^V 2 ] (H(z,z,y)-1) 

= H(p,p,y) (C.l) 

As already noticed in section 4, in eq.( p.l| ) the two dimensional momentum vector p was sub- 
stituted by a complex number p whose norm we name \x: 

p = p 4 + ip 5 ; = ^ ( c - 2 ) 

while p denotes the radius in M. 2 . In section || we dealt with the case in which |c?7(z)| 2 depended 
only on p. In this situation we could expect the integral H (z, z, y) to depend only on p and 
r since the source term and the coefficients of its differential equation would depend only on 
these radial coordinates. Since 7(2) is required to be analytic, this occurs when 7(2) is either a 
logarithm or a simple power of z: 

7M~{'°. g(2) (C.3) 
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Defining the Fourier transform of |c?7(z)| 2 : 

dz A dz 



j(p,p) 



(2vr) 2 



exp [i(pz +pz)/2] \dj(z)\' 



(C.4) 



C being the angular variable on Mr. In the two special cases ( C.3 ) the Fourier transformed source 
term j is a function only of p since we can write: 



i 



(2tt)2 
which impliesf ! 

J(M) = 



2tt 



1 



d( exp [ippt cosC] |#7| 0) = — / pdpJ (pp) \dj\ (p) 



2tt 



lim a ^ ^ jo°° M> J o 

h Io° P d P J o M p=xx* 



3 (log) 



-> oo 
r(q) 



7r(/i 2 )° ! r(l-a) 



(C.5) 



(C.6) 



and where we used the following integral representation of the Bessel Jo function: 

r-27r 



/ dQ exp [i x cos £] = 2ir J (x) 
Jo 



(C.7) 



Under the hypothesis ( p.3|) we can denote the unknown function H by -ff (p, r) and its Fourier 
transform H(p, r) will be given by eq. (4.5). In all other cases H = H(p, p, r). 

However, independently from the explicit form of the twisted field j(z), after partial Fourier 
transform the original equation (4.4) reduces always to: 

L r [a,p]H = J(r,p,p) (C.8) 



which generalizes eq. (4.6) to a more general choice of 7(2) and where the differential operator 
L r [a, p] was defined in eq. Q4.7|). In eq. (|C.8|) we have defined the complete source function: 



J(r,6,p) 



2a 2 

•-0-2 i(p, P) +S(p, r) 



(C.9) 



The Wronskian for a basis of two independent solutions u\(a, p, r) and 1*2(0, p, r) of the ho- 
mogeneous equation ( |4.10 ) associated with eq. ( C. 



is: 



W(ui,U2) = det 



tti(r) u 2 (r) 

u 'l( r ) u 2( r ) 



const x 



p(r) 



const x 



r 3 o(r) J 



(CIO) 



where the overall constant can be adjusted by choosing the overall normalization of the solutions 
u\ and U2- 

Focusing on the case in which the source term J depends only on p, r, the non-homogeneous 
equation (|C.8|) can be solved, according to eq. ( 4.11 ), in terms of the radial Green function 
G (a, p\r,r') defined by eq. (4.12) plus the the boundary conditions to be imposed on the 

10 The Fourier-Bessel transform of |97| 2 in the logarithmic case is a divergent constant, independent from /x. 
Formally we name it j(log) to remind ourselves what it is. At the end, one takes the inverse Bessel- Fourier 
transform and j(log) is restored to 1/p 2 
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solutions H. Following the general theory |3jJ the Green function is constructed in terms of a 
basis ui,ii2 of solutions of the homogeneous equation as follows: 

G( r r') = { U< = °i( r ') u i(r) + c 2 (r') u 2 (r) for r < r' _ 
\u > = di(r')ui(r) + d2{r')u2(r) for r > r' 

where the coefficients q, di become functions of the second variable and are determined by: 

1. the continuity of G at r = r'\ 

u< (r')=u > (r') (C.12) 

2. the discontinuity of l/w(r) in the derivative dG at r = r': 

Sa {i G (r ' r/) lr=r '^ ~ i G (ry) lr=r/ -} = (c,13) 

3. The symmetry of the real Green function in the exchange of its arguments: 

G(r,r') = G(r',r) (C.14) 

4. The boundary conditions at the end-points of the interval [a, oo] that we want to impose 
on the function H(r): 

lim H = h a ; lim H = (C.15) 



r-— >oo 



The same conditions have to be imposed on the Green function. 

While the first three items in the above list are mathematical requirements, the fourth encodes 
the definition of the physical system we want to describe. For this reason we have to choose the 
boundary conditions by careful comparison with the flat case which is our only guide at this 
stage. Furthermore, since the solutions u±,U2 will be shown to have singular behaviors at the 
end points of the interval, eq.( C.15| ) has to be understood in the sense of asymptotic regimes. 



As shown in [37], the first two conditions become a linear system of equations for the r'- 
dependent parameters: 

c = c\—d\ ; d = C2 — d 2 (C.16) 

with the general solution: 

/ / \ _ M r ') ./ /n _ m(r') rrir , 

C(r ' {r' fg{r>)W{r') ' {r ' {r'f g{r')W{r') 1 ' 

Recalling eq.( p.lCj ) if we normalize u\ and U2 in such a way that the constant appearing in the 
Wronskian is fixed to 2: 

W( Ul ,u 2 ) = 3 2 . (C.18) 
r 9v) 

then the Green function takes the form Wfl : 



G(r,r') =G s (r,r')+d 1 (r')n 1 (r) + d 2 (r , )n 2 (r) (C.19) 

where: 



^ (u2(r') u\(r) — u\(r') U2(r)) for r < r' 



G s (r,r') = \ , (C.20) 

tor r > r 
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To fix the last bit of information, namely d\ and d 2 , we have to impose a choice of boundary 
conditions. The analysis of sections and allows to determine the required near cycle 
behavior of the integral H functions and consequently both the expressions for the source term 
S(p, r) ( [4.32 ) and for the two functions d\ and d 2 - 



di(r') = -i« 2 (r') ; d 2 {r') = 



(C.21) 



which determine in turn the form of the Green function, given in eq. (4.31). 

In this final part of the appendix we derive some of the main results of the analysis carried 
out in section O], besides describing some general feature of eq. (4.33). 

Equation ( 4.33|) is characterized by two regular singularities of the homogeneous part at the 
points w = 0, 1. The regular singularity w = corresponds to the position of the exceptional 
cycle r = a while the irregular singularity at w = — oo corresponds to the large distance limit 
r — * oo. The homogeneous equation associated with (4.33) can be framed within the class of 



equations known as confluent Heun differential equations. A confluent Heun equation may be 
written in its standard form as follows: 



dw 2 



+ 



7 & 
4p+ - + 

w w — 



d f (Apaw — 
dw w (w — 1 



<7 



f 







(C.22) 



and in general it has two kind of solutions: an angular solution f^ a \p, a, 7, 5, a; z) and a radial 
solution f( r \p, a, 7, S, a; z), defined as solutions of suitable boundary value problems on the 
intervals ]0, 1[ and ] — 00, 0[ respectively. The homogeneous part of eq. ( 4.33 ) is obtained from 
eq. ( p.22[ ) for 7 = 5 = 1, a = k/(2p) + 1, a = k + 2p in the limit p -> 0. Although the 



a = K/{zp) + 1, a 

literature on this confluent limit of the Heun equation is rather limited, for our purposes we 



could restrict our analysis of the solution of (433) to the neighborhood of the two physically 
relevant singular points, namely w = and w = —00 (i.e. to a region close to the exceptional 
cycle and to radial infinity) and discuss the restrictions imposed on the solution in these two 
regions by the boundary conditions on the interval ] — 00, 0[. As far as the singularity w = is 
concerned, in its neighborhood we can write a formal power series expansion for an integral of 
the homogeneous equation: 



H(n, w) 



w 



n=0 



(C.23) 



the determination of the coefficients c n (n, a) through recurrence relations is a standard pro- 
cedure in dealing with differential equations with regular singularities. Substituting the above 
series in the homogeneous equation we obtain: 



L W H = y~] (cn(fi, a) f (a + n) + Cn-x{^, o) fijp, a + n-1) 



n=0 



+ c {n, cr)/ n (/i, a)) w n = 



(C.24) 



where the coefficients fo(cr + m) and / ra >o(A i ) a + m) are determined by the following expansion: 

00 

(a + m)(a + m-l) + (a + m)P 1 (w) + P 2 (fi,w) = f (a + m) + ^ /„(//, a + 



TO W 



n=l 
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/o(<7 + m) 
fl(H, a + m) 
fn>liP, a + m) 



(cr + m) 2 
k — a — m 
— (k + o + m) 



(C.25) 



Hence equation ( p. 24 ) gives recursion relations for the unknown coefficients c n (/i, a). The first 
of them reads cq([i, a) fo(cr) = and yields a non trivial solution only if the indicial equation is 
fulfilled, namely if /o(cr) = a 2 = 0. The coefficient co(//, a) is thus a free parameter in terms of 
which all the other coefficients are expressed. The solution to a generic order n of the recursion 
relations is represented in eq. ( 4.36| ) that we repeat here in a more explicit form: 



Fn{y>, a) 



rco(M) 



F n (pL, cr) 



F n ip, O") 



flip,, a + n-1) 
foia + n-1) 




f 2 ip, a + n 
flip, a + n 
fo{a + n- 





2) 
2) 



. f n -l(p, a + 1) 
. / n -2(A i , CT + 1) 

• UsiP, a + 1) 



Ma + 1) 



/nG", Cr) 

/ n -lO, cr) 
/n-aO, 0") 

cr) 



(C.26) 



Since <r = is a double solution of the indicial equation, from standard theory of differential 
equations it follows that the homogeneous equation will have a regular solution H^ R \p, r) 
represented by a power series expansion around w = whose coefficients are c^(/i) = c n (fi, a = 
0), and an irregular one H^\/j,, r), both these independent solutions are defined in eq. (|4.35|) . 

As far as the particular solution H^ n ~ h \u, r) of the non-homogeneous equation is concerned, 
we can write a formal power series expansion as in eq. flOp . Expanding the source term on 
the right hand side of eq. ( 4.33 ) as in eq. ( 4.4l| ) and substituting both eq.s ( 4.38 ) and ([4.41 ) in 
eq. (|4.33j ) we obtain: 



h{a) w 



4 n -"Wo(o) + 



(1 - w)(l - 2wf 

oo 

(c^ip) fo(n) + c^Hp) flip, n - 1) + . . . 

n=l 

+ c<f- h) (v)f n (fi, 0) - w n = 



(C.27) 



^From the above equation we may deduce the recursion relations for c« n ^ (u) which are solved 
in eq. ( [4.39| ). Such a solution is repeated here in more explicit form: 



ip) 



(-)" 



Fniti 



Fnip) 



uu foti) 

fxin,n-l) / 2 (jU, n 

/o(ra-l) hip, n 
Mn- 



Fnjp) 

UU r 

2) ... 
2) ... 











fn-lifJ,, 1) 
/n-20, 1) 
fn-sifi, 1) 

/o(l) 



-hifx) s n -i 
-hifj.) S n -2 

-KP) s n-3 



(C.28) 
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Finally let us consider the limit \i S> 1/a, <C 1 and < 1 in which the above series 

expansions can be summed order by order in 1/k according to the expansions ( 4,46j ). In eq.s 
( [4.47D we gave the sum of the series at the first non trivial order. 

For completeness we give below the sums at the next order in 1/k: 



H[ n ~ h \t) 



£ being \kw\. 



R ^ (-2 + n) (-n + n 2 ) f 
c o 2^ 3 ^12 

n=0 



Jt 



£3/2 



^ [2/0(2 ve) (-i+£M£))+\/e (-(Ji(2Ve) (6 + (2+0 io g (o)) + 

£ (2 if 3 (2 v1) + 1 3 (2 >/£) log(0' 
-3 + 2n-3n 2 + n 3 



E 

ra=0 

%) 



so 



Si 



ra! 2 



24 + 24 e - 24 / (2 - 2 ^ I (2 7e) 



+2 V / e/ 1 (2 v / e) + e f /i(2v^) 



(C.29) 
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